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Energy Spectra of the Soft X-ray Diffuse Emission in Fourteen Fields 

Observed with Suzaku 
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Abstract 

The soft diffuse X-ray emission of twelve fields observed with Suzaku are presented together with two 
additional fields from previous analyses. All have galactic longitudes 65° < i < 295° to avoid contributions 
from the very bright diffuse source that extends at least 30° from the Galactic center. The surface bright- 
nesses of the Suzaku nine fields for which apparently uncontaminated ROSAT All Sky Survey (RASS) 
were available were statistically consistent with the RASS values, with an upper limit for differences of 

in i?^5-band. The Ovil and Ovin intensities are well correlated to each other, 



17 x f0~ 6 c s" 1 amin -2 



and Ovn emission shows an intensity floor at ~ 2 photons s 



1 cm- 2 



str -1 (LU). The high-latitude Ovin 



emission shows a tight correlation with excess of Ovn emission above the floor, with (Ovin intensity) = 
0.5 x [(Ovil intensity) - 2 LU], suggesting that temperatures averaged over different line-of-sight show 
a narrow distribution around ~ 0.2 keV. We consider that the offset intensity of Ovn arises from the 
Heliospheric solar wind charge exchange and perhaps from the local hot bubble, and that the excess Ovn 
(2-7 LU) is emission from more distant parts of the Galaxy. The total bolometric luminosity of this galactic 
emission is estimated to be 4 x 10 39 erg s _1 , and its characteristic temperature may be related to the virial 
temperature of the Galaxy. 

Key words: Galaxy: disk — Galaxy: halo — X-rays: diffuse background — X-rays: ISM 



1. Introduction 

The soft X-ray diffuse backgound below I keV is con- 
sidered to consist of the contribution of faint extragalac- 
tic sources and emissions from highly ionized ions, such 
as Cvi, Ovil, Ovin, Fexvn, and Neix, in solar neighbor- 
hoods and in our Galaxy. The extragalactic contribution, 
which we refer to as the Cosmic X-ray Background (CXB) 
in this paper, is estimated to be about 40 % of emission 
in the so-called ROSAT R45 band (McCammon et al. 
2002), which is approximately ~ 0.44— I keV. The emis- 
sion from highly ionized ions is considered to arise from 
at least three different origins. Among them, the solar 
wind charge exchange (SWCX) induced emission from the 
Hcliosphere (Cox 1998; Cravens 2000; Lallement 2004), 
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and the thermal emission from the hot gas in local hot 
bubble (LHB) (McCammon & Sanders 1990) are hard 
to be separated from each other using emission spectra 
above ~ 0.4 keV. In this energy range, the sum of the two 
emission components are approximated by a thin-thermal 
emission of kT ~ 0.1 keV without absorption (Smith ct al. 
2007; Henley et al. 2007; Galeazzi et al. 2007; Kuntz & 
Snowden 2008; Masui et al. 2009). The remaining emis- 
sion is considered to arise from more distant part of the 
galaxy; mostly above or beyond the bulk of absorption 
in the Galactic disk. Kuntz & Snowden (2000) called 
this component the "transabsorption" emission (TAE) 
and separated it in the ROSAT all sky map utilizing the 
directional dependence of the absorption column density. 
They found that the emission spectrum can be described 
by a two-temperature thermal emission model of temper- 
atures kT ~ 0.10 and 0.25 keV. However, because there 
is no constraint on distance other than absorption, it is 
hard to constrain the origins conclusively. 

New insight has been obtained from combined analyses 
of the absorption lines observed in the energy spectra of 
extragalactic objects and emission lines of the same ion 
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species observed nearby. Yao et al. (2009) analyzed the 
absorption spectra of LMC X-3 obtained with the trans- 
mission grating (HETG) onboard Chandra, and the emis- 
sion spectra from the blank fields about 30' away from 
LMC X-3 observed with the CCD camera (XIS) onboard 
Suzaku. The joint spectral fit of the data shows the hot 
gas attributed to the TAE component cannot be isother- 
mal. Instead, a thick Galactic hot gaseous disk whose 
temperature and density decrease exponentially from the 
Galactic midplanc can consistently explain the observa- 
tions. They obtained scale heights of 1.4 and 2.8 kpc, and 
the midplane values of 0.31 keV and 1.4 x 10~ 3 cm -3 for 
the temperature and the density, respectively. 

A similar result was obtained for the sight line toward 
Mrk 421, although the the X-ray emission is mostly based 
on the ROSAT All Sky Survey data, thus the emission 
lines are not resolved spectroscopically (Yao & Wang 
2007). 

In this paper, we analyze energy spectra of the soft X- 
ray diffuse emission of twelve different fields observed with 
Suzaku (Mitsuda et al. 2007) and combine them with 
Suzaku results of two other fields in literature. By virtue 
of the good energy response function of the CCD camera 
(Koyama et al. 2007) and the high sensitivity when com- 
bined with the X-ray mirror (Serlemitsos et al. 2007), 
we clearly detected Ovn emission from all the fields, and 
Ovni emission from most of them. 

In this paper we concentrate on the XIS1 data which 
has much larger effective area below 1 keV than other 
XIS sensors because it employs backside illuminated CCD. 
Throughout this paper, we quote single parameter errors 
at the 90 % confidence level unless otherwise specified. 

2. Observations and data reduction 

2.1. Observations and standard data screening 

In Table 1, we show the observations we used in this 
paper. The analysis results of the seven data sets, 
Off Filament (Off-FIL), On Filament (On-FIL), North 
Ecliptic Pole 1 (NEP1), MBM 12 off cloud (M12off), LMC 
X-3 Vicinity (LX-3), MBM 12 on cloud (M12on), and 
Midplane 235 (MP235) have been already published. 

The data reduction done by Yao et al. (2009) for 
LX-3 is consistent with that shown in this section. Since 
Fujimoto et al. (2007) (NEP1) and Smith et al. (2007) 
(M12off) used version 0.7 processed data, we re-analyzed 
the data from the data reduction. The second observation 
of the North Ecliptic Pole direction (NEP2) was made 
about a half year after the first NEP observation. The 
aim points of the two observations arc identical, although 
the roll angle was rotated by about 180°. 

The two fields at the bottom of Table 1, M12on and 
MP235, are special directions compared to other twelve 
fields. M12on is towards the molecular cloud MBM12 
located at ~ 100 pc with a bright Cataclysmic Variable 
in the field of view behind the cloud, and MP235 is 
in the galactic plane at b = 235° with a large (TVh = 
9 x 10 21 cm -2 ) galactic absorption. They therefore require 
spectral models different from those of the other fields in 



this paper. The spectral results from these two directions 
have been already published (Smith et al. 2007; Masui 
et al. 2009) and were analyzed in the same manner as 
used in this paper, including the geocoronal SWCX re- 
moval. We therefore simply adopt their analysis, although 
we have updated their fits to M12on with more recent cal- 
ibration data. For these two fields we will simply show the 
results in tables and figures and will not show the details 
of analysis in the text. 

In all the observations, the XIS was set to the normal 
clocking mode and the data format was either 3 x 3 or 
5x5. The Spaccd-raw Charge Injection (SCI) was on 
for Low latitude 86-21 (LL21) and Low latitude 97+10 
(LL10). We used version 2.0 processed Suzaku data. We 
first cleaned the data using the selection criteria: elevation 
from sunlit/dark earth rim > 20/5 deg, cut off rigidity > 
8 GV. We checked the dependency of the 0.4 - 0.7 keV 
counting rate on the Oxygen column density of the sunlit 
atmosphere in the line of sight using the MSIS atmosphere 
model (see Fujimoto et al. 2007; Smith et al. 2007; Miller 
et al. 2008). We found that the counting rate was con- 
stant as a function of the column density for the cleaned 
data. Thus there is no significant neutral O emission from 
Earth atmosphere in the filtered data. 

2.2. Removal of Point sources 

We then constructed an X-ray image in 0.3 to 2 keV en- 
ergy range. We detected point sources in the 17/8 x 17/8 
field of view for all the data sets except for LL10. We 
removed a circular region centered at the source position 
from the further analysis. The radius of the circular region 
was determined so that the counts from the point source 
outside the circular region becomes less than 3 % of the 
diffuse X-ray emission in 0.3 to 1 keV energy range. The 
sum of contribution of point sources is estimated to be less 
than 6% of the diffuse emission for all the fields except for 
LH-1 and LL21. The contributions of the two fields, LH-1 
and LL21, are respectively 8 and 12 %. For these fields, 
we analyzed the point source spectra and used the mirror 
point spread function to determine their resiual contribu- 
tion in Ovn and Ovni. We extracted sum of point-source 
spectra from the circular source regions and performed 
model fittings. As the background spectra, we used the 
spectra extracted from source-free circular regions whose 
distances from the optical axis of the telescope are equiva- 
lent to the source regions, respectively. We estimated the 
upper limits of Ovn and Ovni emission intensities of the 
source regions with spectral fits, then their contamination 
to the SXDB spectra. The upper limit of contamination 
was 0.4 LU for Ovn emission for the both fields. It was 
0.2 LU and 0.6 LU for Ovni emission for LH-1 and LL10, 
respectively. For GB1428+4217 (GB) a point source of an 
intensity of 2.6 x 10 -13 erg s" 1 cm -2 in the energy band 
of 0.3 to 1 keV ( 22 x 10~ 13 erg s" 1 cm~ 2 in 1 to 10 keV) 
was detected at the center of the field of view. We re- 
moved a circular region of a 5 arc minute radius for this 
observation. The contribution of the source in the count- 
ing rate of the remaining region is estimated to be only 
0.2 %. 
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Table 1. Log of observations, ordered by |6 



Data set 


Obs ID 


Date 


Expo: 


mre (ks) 


Aim point 


ID 


Field Name (Short Name) 






Total 


Cleaned 


(1,6) 


(-ELon, -^Lat)* 


1 


GB1428+4217 (GB) 


701092010 


Jun 12-13, 2006 


48.7 


34.9 


(75.9, 64.9) 


(194.2, 52.7) 


2 


High latitude B (HL-B) 


500027020 


Feb 17-20, 2006 


103.6 


29.7 


(272.4, -58.3) 


(4.4, -61.4) 


3 


Lockman hole 2 (LH-2) 


101002010 


May 17-19, 2006 


80.4 


40.0 


(149.7, 53.2) 


(137.1, 45.1) 


1 


Lockman hole 1 (LH-1) 


100046010 


Nov 14-15, 2005 


77.0 


61.7 


(149.0, 53.2) 


(137.2, 45.5) 


5 


Off Filament a (Off-FIL) 


501001010 


Mar 1-2, 2006 


80.1 


59.6 


(278.7, -47.1) 


(354.8, -72.6) 


6 


On b llament (On-i 1 1L) 


501002010 


Mar 3-6 , 2006 


101.4 


59.2 


(278.7 , -45.3) 


(354.1, -74.4) 


7 


High latitude A (HL-A) 


500027010 


beb 14-15, 2006 


73.6 


53.2 


(68.4, 44.4) 


(228.8, 63.5) 


8 


MBM12 off cloud b (M12off ) 


501104010 


Feb 6-8, 2006 


75.3 


51.0 


(157.3, -36.8) 


(44.5, 2.3) 


9 


LMC X-3 Vicinity c (LX-3) 


500031010 


Mar 17-18, 2006 


82.0 


56.1 


(273.4, -32.6) 


(41.2, -86.2) 


10 


North Ecliptic Pole l d (NEP1) 


100018010 


Sep 2-4, 2005 


106.2 


58.7 


(95.8, 28.7) 


(334.8, 88.7) 


11 


North Ecliptic Pole 2 (NEP2) 


500026010 


Feb 10-12, 2006 


75.6 


16.5 


(95.8, 28.7) 


(334.8, 88.7) 


12 


Low latitude 86-21 (LL21 ) 


502047010 


May 9-10, 2007 


81.5 


57.0 


(86.0, -20.8) 


(347.6, 38.4) 


13 


Low latitude 97+10 (LL10) 


503075010 


Apr 15-16, 2008 


79.8 


40.8 


(96.6, 10.4) 


(0.7, 70.6) 


Rl 


MBM12 on cloud b - c (M12on) 


500015010 


Feb 3-6, 2006 


102.9 


68.0 


(159.2, -34.5) 


(47.2, 2.6) 


R2 


Midplane 235 c (MP235) 


502021010 


Apr 22-25, 2007 


189.5 


53.0 


(235.0, 0.0) 


(119.5, -40.6) 



Results previously published by a Henley et al. (2007), b Smith ct al. (2007), c Yao et al. (2009), d Fujimoto et al. 
(2007), c Masui et al. (2009). 
* Ecliptic coordinate 



2.3. Removal of geocoronal SWCX 

The spectrum beiow 1 keV could be contaminated by 
the SWCX induced emission from the geocorona (e.g. 
Cravens et al. (2001)). Fujimoto et al. (2007) reported 
a flarc-likc increase of X-ray flux during the Suzaku ob- 
servation of NEP1, which was clearly visible in the light 
curve of 256 s time bins. In all other observations, the 
increase is too small to be recognized in the light curve of 
256 or 512 s time bins, but still it could be statistically 
significant in the X-ray spectra when they are integrated 
over a certain amount of time (e.g. > 30 ks). Thus as the 
last stage of the data reduction, we removed time inter- 
vals in which the data can be contaminated by the SWCX 
from the geocorona. 

We consider two parameters which are related to 
Suzaku's geocoronal SWCX inensity. The first one is 
the solar wind flux near the Earth. We calculated the 
solar- wind proton flux at the Earth using ACE SWEPAM 
data 1 . When the ACE SWEPAM data are not available, 
we used WIND SWE data 2 or OMNI data from CDAWeb 
(Coordinated Data Analysis Web) 3 . 

The second parameter is the Earth-center to magne- 
topause (ETM) distance, where the magnetopause is de- 
fined as the lowest position along the sight line of Suzaku 
where geomagnetic field is open to interplanetary space. 
The ETM distance is in the range of ~ 2 to ~ 10 Earth 
radii (Re ), while Suzaku is in a low earth orbit of an 
about 650 km altitude. Thus it is always looking at the 
sky through the magnetopause. The probability of a con- 
tamination in Suzaku spectra by the SWCX from the 
geocorona increases if the shortest Earth-center to mag- 
netopause (ETM) distance is < 5 Re (Fujimoto et al. 
2007). We calculated the ETM distance every 256 s for 
all the observation periods using the T96 magnetic field 



1 The data available at http://www.srl.caltech.edu/ACE/ACS/ 

2 The data available at 
http: / /web. mit.edu/afs/athena/org/s/space/www/wind. html 

3 The data available at http:/ /cdaweb. gsfc.nasa.gov/cdawcb/sp_phys 



model (Tsyganenko & Sitnov 2005) . We obtained the 
interplanetary plasma parameters required for the calcu- 
lations from the CDAWeb. 

When the sun is quiet the solar wind proton flux stays 
in the range of (1 — 4) x 10 s protons s _1 cm -2 , and shows 
slow (times scales of ~ 10 ks) time variations. On the 
other hand, during flaring events, it increases and can go 
up to ~ 5 x 10 9 protons s _1 cm -2 , and shows fast (< 1 ks) 
time variations. 

For the four observations, GB1428+4217 (GB), 
Lockman hole 1 (LH-1), On Filament (On-FIL) and Low 
latitude 86-21 (LL21), the proton flux was always be- 
low 4 x 10 8 protons s -1 cm~ 2 . Among them, during 
GB and On-FIL observations, the proton flux showed 
slow variation and the data sets can be divided into 
two subsets of (1 — 2) x 10 9 protons s _1 cm -2 and (2 — 
3) x 10 9 protons s' 1 cm~ 2 for GB, and of (2 - 3) x 
10 9 protons s _1 cm~ 2 and (3 — 4) x 10 9 protons s _1 cm' 2 
for On-FIL, respectively Both during the GB and On- 
FIL observations, the ETM distance varied in 2 to 10 Re 
range. We created spectra for the subsets and performed 
spectral fits to determine the Ovu and Ovni emission in- 
tensities. The fitting model used in this analysis will be 
described in detail in the next section (3.1). We found 
the difference in the Ovil line intensities was less than 2 
cr for both fields. This corresponds to 1.7 LU and 1.3 LU 
for GB and On-FIL, respectively. For Ovni, the differ- 
ence was within 1 a, which is 0.5 and 0.6 LU for GB and 
On-FIL, respectively. 

Among all other observations, in part of HL-B, HL- 
B, and LH-2 observations, the proton flux stayed be- 
low 4 x 10 8 protons s _1 cm~ 2 , and showed variations of 
more than 1 x 10 8 protons s _1 cm -2 was observed. We 
thus also subdivided the data sets according to the pro- 
ton flux. We found that the variations of Ovu and 
Ovni intensities were within 2 a level for all those data 
sets. We thus conclude that when the proton flux is 
<4x 10 8 protons s _1 cm~ 2 the contamination of geocoro- 
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nal SWCX is relataively small, although there still could 
remain possibility of contamination of ~ 1.5 LU and ~ 0.5 
LU levels for Ovu and Ovili, respectively. We thus de- 
cided to use all the data for the four observations, (GB, 
LH-1, On-FIL), and LL21, 

For the seven observations, High latitude B (HL-B), 
Lockman hole 2 (LH-2), Off Filament (Off-FIL), LMC X- 

3 Vicinity (LX-3), North Ecliptic Pole (NEP1), the second 
North Ecliptic Pole observation (NEP2) and Low latitude 
97+10 (LL10), the proton flux shows shows flarc-likc in- 
crease and exceeds 4 x 10 8 protons s _1 cm -2 in about a 
half of the obsevation. We thus subdivided those data 
into two subsets according to the proton flux, and deter- 
mined the Ovu and Ovin intensities. We found signifi- 
cant (more than 2 a in Ovu or Ovin) difference between 
the two spectra for four cases: HL-B, LH-2, NEP1 and 
NEP2. An example of spectral comparisons is shown in 
Figure 1. For these four data sets, we decided to use only 
the time intervals in which the proton flux was lower than 

4 x 10 8 protons s _1 cm~ 2 . These results are consistent 
with the geocoronal SWCX contamination reported for 
Suzaku observations of 4U1830-303 vicinities (Mitsuda et 
al. 2008). 

In order to find out why no difference was found in the 
remaining three observations, we checked the ETM dis- 
tance of these observations. As already reported in Yao 
et al. (2009), the ETM distance became as short as 1.4i?E 
during the LX-3 observation. However, the magnetic field 
was open to anti-sun direction when the ETM distance 
was < 5i?E during which solar-wind particles cannot pen- 
etrate. This explains why there was no significant differ- 
ence between the high and low proton-flux spectra for this 
data set. We thus decided to adopt all the data for this 
observation. During the LL10 and Off-FIL observation, 
the magnetic field was occasionally open to sun side when 
the ETM distance dropped below 5i?E- We decided to 
simply discard the time intervals with high proton fluxes 
for these data because we still have enough statistics. 

There remain two data sets, High latitude A (HL- 
A) and MBM12 off cloud (M12off). During the HL- 
A observation, the proton flux varied in the range of 
3 — 6 x 10 8 protons s _1 cm~ 2 but the time interval ex- 
ceeding 4 x 10 s protons s _1 cm~ 2 was short (~ 1/8 of 
data). During the M12off observation, it stayed at a level 
of ~ 5 x 10 8 protons s _1 cm~ 2 . Thus for these data sets, 
we can not compare spectra with low and high proton 
fluxes. For HL-A, we subdivided the data sets accord- 
ing to the ETM distance and created two energy spectra 
with ETM distances = 5 - 10i? E and > 10R E . We found 
that the two spectra with different ETM distances show 
no significant difference. Thus we will use data with the 
ETM distances > 5i?E- There was a small fraction of the 
data with the ETM distance < 5 -Re, which we decided to 
discard. For M12off, the ETM distance stayed in 7- 9i? E - 
Thus further subdivision was not possible, and we decided 
to use all the data. 

We tried to exclude time intervals for which contami- 
nation of time variable geocoronal SWCX was significant. 
However we should keep it in mind that there is still un- 




1 2 
Energy (keV) 



Energy spectra of High Latitude B (HL-B) 
time intervals with solar wind proton flux 
squares) or lower (open squares) than 
ons s — 1 cm -2 . The vertical error bars arc the 



Fig. 1. 

for the 
higher (filled 
4 X 10 s proti 

1 — a statistical errors. The non X-ray background was sub- 
tracted. There are large (> 3 <r) discrepancies between the 
Ovu emission (~ 0.56 keV) intensities of two spectra, which 
indicates significant geocoronal SWCX contamination in the 
higher proton-flux data. 

certainty of contamination by geocoronal SWCX at 1.5 
LU and 0.5 LU levels for Ovu and Ovin, respectively. 
In addition, the last data set, M12off, could have larger 
contamination. In Table 2, we summarize the process we 
applied to remove time intervals which are suspected to 
be contaminated by the SWCX from the geocorona. We 
also show the process for M12on and MP235 fields in the 
table. 

2. 4- Non X-ray background 

The non X-ray background (NXB) spectrum was con- 
structed from the dark Earth database using the standard 
method in which the cut off rigidity distributions of the 
on-source and the background data were made identical 
(Tawa et al. 2008). We found 5 to 10 % discrepancies 
in the counting rates above 10 keV between the NXB and 
the observation data, suggesting background uncertainty 
of this level. We also analyzed front-iluminated CCD data 
(XIS0, XIS2, and XIS3. but XIS2 was not available af- 
ter November 9, 2006) to determine the intensity of the 
Cosmic X-ray Background (CXB) above 2 keV for each 
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Table 2. Summary of Geocoronal SWCX removal process 



ID 


(1) 


(2) (3) (4) (5) (6) 


1 (GB) 


/ 




2 (HL-B) 




✓ 


3 (LH-2) 




✓ 


4 (LH-1) 


✓ 




5 (Off-FIL) 




/ 


6 (On-FIL) 






7 (HL-A) 




/ 


8 (M12off ) 




/ 


9 (LX-3) 




/ 


10 (NEP1) 




✓ 


11 (NEP2) 




✓ 


12 (LL21 ) 


✓ 




13 (LL10) 




/ 


Rl (M12on) 




/ 


R2 (MP235) 




/ 


(1) The solar wind flux was 


alway < 4 x 10 s protons s _i cm -2 . All time intervals were 


adopted. 
(2) Excess in 


the high-solar-wind-flux spectrum. Time intervals with high solar wind flux 



were discarded. 

(3) The geomagnetic field was open to the anti-sun side, when the ETM distance was low 
(< 5Re)- All time intervals were adopted. 

(4) The geomagnetic field was occasionally open to the sun side when the ETM distance 
was low. Time intervals with high solar wind flux were discarded. 

(5) Proton flux stayed in 3 — 6 x 10 s protons s" 1 cm~ 2 . Time intervals with ETM distance < 
5i?E were discarded. 

(6) Proton flux stayed in 3 — 7 x 10 8 protons s _1 cm" 2 . ETM distance was always > 10i?E- 
All time intervals were used. 



sensors. We determined the CXB intensity by the spec- 
tral fits. The model function will be described in section 
3. The intensities determined by different sensors were 
consistent with each other within the 90% statistical er- 
rors for all observations except for NEP2. For this obser- 
vation, the intensity from XIS1 were not consistent with 
those from other three sensors. It was necessary to adjust 
the NXB level by about 30% for XIS1 to obtain consis- 
tent result. This suggests that not only the NXB level but 
also its spectral shape may vary from the one extracted 
from the standard NXB data base. The CXB intensity 
of NEP2 after the NXB correction was consistent with 
the CXB intensity estimated from NEP1 data. Since in 
the energy range below 1 keV, the non X-ray background 
is only about 10 % of the diffuse X-ray background, the 
spectral parameters derived in section 3 are not sensitive 
to this level of background uncertainty. 

In Suzaku NXB non-subracted energy spectra, there are 
three instrumental lines at 1.486, 1.740, and 2.123 keV 
(Tawa ct al. 2008). Since the line intensity and the line 
spread function of the detector vary with time, the resid- 
uals of those lines appears in the NXB subtracted spectra 
as lines or line-like structures in the 1.3 to 2.1 keV range. 
In order to study those features in our spectra, we first 
fitted all the NXB subtracted spectra in the energy range 



1.1 to 5 keV. We employed a power-law function absorbed 
with a Galactic absorption for the continuum and narrow 
Gaussian functions for the NXB lines. We found that 
we need to include a line at 1.828 keV for HL-B, and a 
line at 2.157 keV for LL10 to obtain an acceptable fit. 
For other data sets, such a line was not required. In the 
further spectral fits, we include these Gaussian functions 
with the parameters fixed at the best fit values for these 
two data sets. The results of the further analysis do not 
depend on inclusion or non-inclusion of these lines in the 
models of spectral fits. 

2.5. Arf and rmf files 

In order to perform spectral fitting, we generated 
an efficiency file (arf file) for a spatially uniform 
emission, using the xissimarfgen software version 2008- 
04-05 (Ishisaki et al. 2007), assuming a 20'-radius 
flat field as the input emission of the generator. A 
pulse height re-distribution matrix (rmf file) was 
created by the script xisrmfgcn version 2007-05-14. 
The calibration files, ae_x;il_quanteff_20070402.fits, 
ae_xil_rmfparam_20080311.fits, 

ae_xil_makepi_20080131.fits, and 
ae_xil_contami_20071224.fits were used to generate 
the arf and rmf files. 
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2. 6. Contamination of the optical blocking filter 

The degradation of low energy efficiency due to the con- 
tamination on the XIS optical blocking filter was included 
in the arf file. The gradient of the contaminant thick- 
ness over the optical blocking filter is also taken into ac- 
count. Systematic errors in the contaminant thickness 
are estimated to be about 10 % (The Suzaku Technical 
Description 4 , see also Fujimoto et al. (2007) for early 
data and Yamasaki et al. (2009) for recent data). We 
performed all the analysis described in the next section 
not only for the nominal contaminant thicknesses but also 
assuming 1 0% thicker and thinner contaminants. Because 
the contaminant thickness is gradually increasing with 
time, we created the arf files of 10% thicker or thiner con- 
taminants by shifting the observation dates from the real 
dates, except for the 10 % thicker cases of LL21 and LL10. 
For the two cases, the +10 % thickness is larger than that 
of the most recent date in the latest calibration database. 
Wc thus added an absorption model in the model spectra 
to represent the extra contamination. For all data sets, 
we found that the best-fit parameters changed by only 
small amount when we varied the assumed contamination 
thickness by ±10 %. Among the parameters, the Ovn 
emission intensity is most sensitive to the contaminant 
thickness. The Ovn emission intensity varied at most by 
0.8 LU (LU — photons s _1 cm~ 2 str~ 1 ), while the statisti- 
cal errors were 0.5 to 1. 5 LU. Thus we will show only the 
results with nominal contamination thickness. 

3. Analysis and results 

3.1. Determination of Ovn and Ovni emission intensi- 
ties 

In Figure 2, the energy spectra of thirteen data sets 
are shown, together with the reference data set, M12on. 
(The spectrum of MP235 is shown in Figure 5). We can 
clearly see the Ovn (~ 0.56 keV) emission lines in all of 
the spectra. 

We first fitted the spectra in 0.4-5.0 keV with a model 
function, consisting of the CXB component represented 
by a power-law function or broken power-law functions 
absorbed with Galactic neutral medium, and two thin 
thermal-emission components which we consider to rep- 
resent the Heliospheric SWCX induced emission plus the 
thermal emission from the local hot bubble, and the ther- 
mal emission from more distant part of the Galaxy. As the 
thin thermal emission model, we used APEC (http://hea- 
www.harvard.edu/APEC/). The first APEC component 
representing the sum of the SWCX from the Hcliosphere 
and the emission from the local hot bubble (hereafter 
SWCX+LHB) is unaffected by the Galactic absorption. 
The contribution of the LHB is considered to be at most 
one third of this component and could be as small as 
~ 1/20 (see the discussion by Smith et al. (2007)). 
The SWCX induced X-ray emission is a non-thermal pro- 
cess. Thus there is no logical reason that the spectrum 
of this component can be approximated by thermal emis- 

4 available from http:/ /www. astro. isas.jaxa.jp/suzaku/doc/suzaku. 
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sion. In particular, fine structure of lines and line in- 
tensity ratios of different transitions must be very differ- 
ent from those of thermal emission (e.g. Kharchcnko ct 
al. (2003)). However, the spectrum of this component 
has been successfully represented by thin thermal models. 
Since Suzaku does not have enough energy resolution to 
resolve fine structures of lines, and since this component 
mainly contributes to the Suzaku spectra of the energy 
range below ~ 0.6 keV, the same approximation can be 
applied to the present data. 

Galactic emission beyond the LHB may arise inside the 
bulk of Galactic absorption. However, we assume that the 
second thermal component is subject to the total Galactic 
absorption; i.e. the transabsorption emission (TAE). This 
assumption is consistent with the discussion in the next 
subsection (3.2). We fixed the column densities of absorp- 
tion, Ah, for the CXB and the TAE components to the 
total H column density determined by 21cm observations 
(Dickey & Lockman 1990) except for On Filament (On- 
FIL) and Off Filament (Off-FIL) data. For these two data 
sets, Henley et al. (2007) and Henley & Shelton (2008) 
used the column densities estimated from the 100-^m in- 
tensity from the all-sky IRAS maps using the conversion 
relation for the southern Galactic hemisphere by Snowdcn 
et al. (2000). We thus tried spectral fits using two differ- 
ent Ah values. 

The average AGN spectrum below ~ 1 keV becomes 
steeper than above ~ 1 keV, and the average photon in- 
dex was determined to be 1.96 below ~ 1 keV by Hasinger 
et al. (1996). We thus tried two different models for the 
CXB: a power-law model with a photon index fixed to 
1.4, and a double broken power-law model used by Smith 
et al. (2007). The latter model consists of two broken 
power-law functions with photon indices of 1.54 and 1.96 
below 1.2 keV and a photon index of 1.4 above the en- 
ergy. We fixed the normalization of the broken power law 
component with a low-energy photon index of 1.54 to 5.7 
photons s _1 keV -1 str -1 at IkeV and set the normaliza- 
tion of the other broken power-law component free. 

First we fixed the elemental abundances of both thin 
thermal components to the solar value (Anders & Grevesse 
1989), and set the temperature and normalization free. 
We call these models with the double broken power-law 
CXB model and with the simple power-law CXB model 
respectively model 1 and model 1'. 

In Tables 3 and 4, we summerize the results of the spec- 
tral fits. The best fit model functions convolved with the 
instrument response functions are plotted in Figure 2 for 
fits with model 1. For HL-B and HL-A with model 1, and 
for GB and HL-B for model 1', the temperatures of the 
two thermal components could not be constrained well be- 
cause of strong coupling between the two components. In 
these cases, we fixed the temperature of the SWCX+LHB 
components to the average values of other spectra. For 
GB and LH-2 with model 1, the existence of the TAE 
component was not significant. Thus for these two cases 
we fixed the temperature of the TAE component to the av- 
erage value of other spectra other than LL10 (0.264 keV), 
_tdyhd estimated the upper limit of the normalization pa- 
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Fig. 2. Observed spectra (crosses) with NXB subtracted, best-fit model 1 and its components (step functions) convolved with the 
instrument response function and residuals of the fit (bottom panels). The model function consists of three spectral components, 
CXB, TAE and SWCX+LHB. The TAE and SWCX+LHB were respresented by thin thermal emission model of solar abundance. 
The CXB component was represented by a double broken power-law model. Thus there are two curves for this component. The 
vertical error bars of data points correspond to the 1 a statistical errors. Some of the spectra show a large excess over the model at 
Ne K emission (~ 0.9 keV). The hard spectral component which is dominant above ~ 2 keV for M12on represents the Cataclysmic 
Variable behind the MBM-12 molecular cloud. The observed spectrum of MP235 is found in Figure 5. 
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Fig. 2. Continued. 
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Table 3. Results of three-component (CXB, TAE, SWCX+LHB) spectral fits with the double broken power-law CXB (model f) 



CXB 5 TAE SWCX+LHB X 2 /dof 
ID (10 20 cm~ 2 ) Norm c kT (kcV) Norm d kT (keV) Norm d 



1 (GB) 


1.40 


5.8±o.e 


0.264° 


< 1.7° 


1 42 +u - U4J 
u.i4z_ 020 


6.2+1.4 


91.25/92 


2 (HL-B) 


3.36 


2.7±o.6 


n 979 +0.081 
U.Z(Z_q 076 


1 6+ 16 

L -°-0.S 


0.120 f 




94.45/87 


3 (LH-2) 


0.56 


3.6±o.5 


0.264° 


< 1.0° 


0.099 +0.027 


16.31^9 


94.74/99 


4 (LH-1) 


0.56 


5 fi+ - 4 

°- D -0.5 


47fi+ - 279 

u -*'°-0.149 


1 0+0.8 

L - 6 -0.5 


n ,07+0.011 

U.IO ( _o.017 


8.3I 3 - 7 


154.08/127 


5 (Off-FIL) 


4.19 


3.0±o.4 


n 971 +0.018 
U.Zf 1_ 023 


' - z -i.o 


n 1 18 +0.006 
U - 11S -0.008 


25.5+^ 


173.17/122 


5' (Off-FIL) 


1.90 


2.9±o.4 


PI 97Q+0.011 
U.Z/o_ 024 


5 8+ 1 - 4 

°-°-0.9 


n 199+O.OI5 

u.izz_ 009 


21.9±lj 


169.90/122 


6 (On-FIL) 


4.61 


4.8±o.5 


R44+ - 081 

U.D++_ o7 5 


1 q+0-5 


0.182±S;8« 


7.0+0.7 


128.57/117 


6' (On-FIL) 


9.60 


c 1+0.5 
°' 1 -0.6 


627+ 072 

u - oz ' -0.063 


9 7+O.6 
z -'-0.5 


0.183±o.oio 


1+0.6 

X- L -0.7 


127.35/117 


7 (HL-A) 


1.02 


4.9+0.5 


n 907+O.045 
U.ZO(_ 039 


3.511;? 


0.120 f 


10. o+^- 2 


141.01/116 


8 (M12off) 


8.74 


1.8+0.5 


245+ 0104 
u.z+o_ 118 




102+ 016 
u.iuz_ 080 


21.4i 21 5 1 


86.27/89 


9 (LX-3) 


4.67 


8.7+0.5 


n 01 O+0.087 
U.OlU_ 029 


6-0^ 


n i4n+ 031 




211.71/196 


10 (NEPI) 


4.40 


o n +0.5 
d.U_ 4 


2 44+ - 023 

U.Z44_ 017 


9 5+ 1 - 7 

y -°-1.9 


n 1 1 0+0. 009 
u.iio_ 0011 


23.8+^ 


158.44/120 


II (NEP2) 


4.40 


2.8+0.8 


n 971+0.040 
u - z ' J^-o.oso 


7-7^3 


q llfi +0.020 
U.llO_ 025 


18.9t 8 7 5 3 


60.85/55 


12 (LL21) 


7.24 


4.3+0.5 


0-284±g;g« 


6-8±?; 7 


109+ 019 
u.iuy_ .020 


24.5t 27 6 ° 


105.03/97 


13 (LL10) 


27.10 


2 6 +0 - 7 

z -°-0.6 


0.738l°; 2 g 


1 c+0.6 

i '°-0.7 


119+0.049 

U.11Z_ .032 


1+26.2 


84.28/ 89 


Rl M12on« 


40.00 


1.7+0.6 






0.095±oT 7 


i8.8±g:8 


141.62/136 


R2 MP235 


90.00 


6.0+1.0 


( 0.765l8;°g 


3.73+0.38 ) h 


o.io8l8:S? 


12.4i 6 7 4 9 


75.80/84 



a The absorption column densities for the CXB and TAE components were fixed to the tabulated values. 
They were estimated from Dickey & Lockman (1990), except for 5' and 6', for which the values were 
estimated from the 100 /xm intensity of the direction (Henley & Shelton 2008). 

b Two broken power-law model was adopted. The photon indices below 1.2 keV were fixed to 1.52 and 
1.96. The normalization of the former index component is fixed to 5.7 phtons s~ cm" " 2 keV" 1 str- 1 @lkeV. 
° The normalization of the broken power-law component with a photon index of 1.96 below 1.2 keV. 

d The emission measure integrated over the line of sight, (l/47r) J n e n\{ds in the unit of 10 14 cm~ 5 str -1 . 
e The TAE component was not significant. The upper limit of the normalization was determined for the 
average TAE temperature of the TAE component of other observations. 

f Because of the strong coupling between the TAE and LHB+SWCX components, the temperatures of the 
two components were not well determined. We thus fixed the temperature of the SWCX+LHB component 
to the average value of other obervations. 

s A spectral component for the bright point source represented by an absorbed brcmsstrahlung was included 
in the fits but is not tabulated. The model parameters are kT — 200 keV(fixed), the normalization factor 
= 2.66 ± 0.14 x 10~ 12 cm~ 5 str -1 (when a point source was replaced with a 20' radius flat emission), and 
A+i = 5.5 x 10 22 cm- 2 (fixed). 

h This component represents the "narrow bump" component (Masui et al. 2009) rather than the TAE. No 
absorption was included in this component. This component emits Ovin emissions and the O abundance 
was set free (the best-fit value is 3.1 1}' 3 solar. 
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Table 4. Results of three-component (CXB, TAE, SWCX+LHB) spectral fits with the power-law CXB model (model 1') 





A H a CXB b 


TAE 


SWCX- 


-LHB 


X Vdof 


ID 


(10 20 cm- 2 ) Norm c 


kT (keV) Norm d 


kT (keV) 


Norm d 





1 (GB) 


1.40 


10.9±o.6 


n (Toe 

u.zuo_ 050 


n — t 

q 0+4.1 

o.o_! 4 


0.112 6 


TTrTj 

8-i-4: 7 


89.86/91 


2 (HL-B) 


3.36 


8.0±o.5 


n 97'} : 0.006 
u - z '°-0.043 


9 q+1.0 
z -°-0.9 


0.112° 


9 5 +2 - 6 


96.37/87 


3 (LH-2) 


0.56 


8-7±g;l 


fl qqc+0.160 

u.ooo_ 148 


1 n +1 - 7 


0.087^;EJ 2 


36.4+^ 8 9 5 


89.76/97 


4 (LH-1) 


0.56 




U.OO*_p 038 


2.5±o.e 


0.1161^08 


15.21 4 ? 


172.90/127 


5 (OFF-FIL) 


4.19 


8.2±o.4 


n 97q+0.017 
u.zfo_ oi9 


7 q+i-5 


U.llO_ oo 6 


28.9±1;| 


178.04/122 


5' (OFF-FIL) 


1.90 


8.2±o.4 


97fi+ 012 
U.Z(D_ 021 


°-°-0.9 


0.118t";oo8 


25.7±g-f 


174.37/122 


6 (ON-FIL) 


4.61 


10.0±g;f 




2.1±0.5 


0.182±o.on 


7-8toi 


124.99/ 117 


6' (ON-FIL) 


9.60 


10.2±o.s 


0-618^;°f 6 


q n+0.5 
°- u -0.6 


n i oq+0.010 
U.lOO_ 009 


Q+0.7 
8 *°-Q.6 


131.00/117 


7 (HL-A) 


1.02 


10.1±0.4 


r\ n a ^.4-0 036 

U.z<±z_ 04g 


a r\-\-H n 

4.9^-0 


u.uao_ 04 6 


O 7 Q ~f~56 . 1 

13 3 


142.17/115 


8 (M12off) 


8.74 


7.2±o.4 




2-7ti.o 


101+ 0015 

U.1U1_ 0X2 


90 q+io:e 


86.99/89 


9 (LX-3) 


4.67 


13.4±o.5 


0.308±g;i 


7 6+ 12 

' -°-1.5 


u.izo_ 013 


14.7±i"f 


214.27/196 


10 (NEP1) 


4.40 


8.1±0.4 


fl 250+ 023 
u.zou_ oi 6 




n -, -, 9+O.OO8 
U - 11Z -0.011 


26.81" 


171.74/120 


11 (NEP2) 


4.40 


8.1±0.7 


0.275±o.o36 


S.St 2 ; 9 


q n o+0.024 
U.llO_ 023 


21.8± 9 9 6 2 


63.15/55 


12 (LL21) 


7.24 


9.5±o.s 


n 987+0.033 


7 6+ 1 - 7 

' - D -1.5 


q lf)1 +0.018 
U.1U1_ 012 


oq r+25.7 


104.43 /97 


13 (LL10) 


27.10 


7 O+0.5 


n 791+0.190 

u - ' zl -0.166 


i q+0-8 


n -, -, n+0.052 
U.llO_ 030 


8-0 + 5 9 3 9 


84.70/89 


Rl M12on f 


40.00 


7 9+0.6 
' - z -0.5 






096 +u ui4 
u.uyo_ 007 


18.5+^ 


141.34/136 


R2 MP235 


90.0 


11.1+0.9 


m 766+ 041 

l u - ,OU -0.039 


q 75+0-40w 

J - '°-0.37i 


U.1UO_ .032 


14 1 +64 


76.54/84 



a The absorption column densities for the CXB and TAE components were fixed to the tabulated values. 
They were estimated from Dickey & Lockman (1990), except for 5' and 6'. The column densities were 
estimated from the 100 fim intensity of the direction. 
b A power-law model with a photon index of 1.4 was adopted. 

c The normalization of the power-law function with the unit of phtons s - cm" " 2 keV _1 str _1 @lkeV. 
d The emission measure integrated over the line of sight, (l/47r) J n c nnds in the unit of 10 14 cm~ 5 str . 
Because of the strong coupling between the TAE and LHB+SWCX components, the temperatures of the 
two components were not well determined. We thus fixed the temperature of the SWCX+LHB component 
to the average values of other obervations. 

f The model parameters for the point source were kT — 200 keV(fixed), the normalization factor = 2.65 ± 
0.14 x 10~ 12 cm~ 5 str -1 (for a 20' radius flat emission), and TVh = 5.5 x 10 22 cm~ 2 (fixed). See table note (g) 
of Table 3. 

s "Narrow bump" component. The best-fit O abundance was 3.1±\i- See table note (h) of Table 3. 



ramctcr. Although some different fitting procedures were 
necessary between models 1 and 1' for a few data sets as 
described above, the best-fit parameter values of the TAE 
and SWCX+LHB components determined using different 
CXB models are consistent within the statistical errors. 
We performed all the fits shown in this paper using two 
different CXB models. However, since all of those results 
are consistent with each other within the statistical er- 
rors, we will only show the results with the double broken 
power-law model for the CXB in the further analysis. 

For the On-FIL and Off-FIL results in Tables 3 and 
4, we find the best-fit parameters to be consistent with 
each other within the statistical errors for both absorption 
values. We performed all the fittings shown later for these 
data sets with both column densities but continued to find 
consistent results, and so will show only the results using 
the absorption column densities from Dickey & Lockman 
(1990). 

The resultant x 2 values were generally good (reduced x 2 
is 0.95 - 1.42 for 87 - 196 degrees of freedom.) The best-fit 
values of the temperatures of the two thin-thermal models 
are contained in a relatively narrow ranges. They are 



from 0.10 to 0.14 keV for the LHB + SWCX component 
and from 0.22 and 0.48 keV for the TAE, except for the 
TAE temperatures of On-FIL (kT = 0.644 keV) and LL10 
(kT = 0.738 keV). 

We then estimated the surface brightness of Ovil K a 
and Ovin K a emission. We set the Oxygen abundances 
of the two APEC models to and added two Gaussian 
functions that represent Ovil K a and Ovin K a emission. 
We set the intrinsic width of the Gaussian functions to a 
value small compared to the detector energy resolution. 
We did not include absorption in the Gaussian compo- 
nents. We set the centroid energies and the normalization 
parameters of the Gaussians and the normalization pa- 
rameters of the CXB and APEC components free. We set 
the intensities of the TAE component to for GB and 
LH-1, because we had obtained only upper limits for this 
component. The surface brightness of a line is calculated 
from the best-fit value of the Gaussian normalization. 

The Ovil K/3 emission is at 666 eV which will there- 
fore blend with the Ovin K Q emission at 654 eV. We 
estimated the Ovil Kp intensity from the best-fit value 
of Gaussian normalization of Ovil K a and subtracted it 
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Table 5. Results of spectral fits to determine Ovu and Ovm line intensities with broken-power-law CXB model (model 1) 
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CXB a 


XAE a,b 


SWCX+LHB b 


Ovu 


Ovm 1 


X Vdof 


ID 


Norm c 


Norm d 


Norm d 

A r+b.6 


Centroid Norm 6 
n r^+u.uiu o P+u.y 


Centroid Norm c 

n or, j+U.UIY r, - I _ . 





2 (HL-B) 

3 (LH-2) 

4 (LH-1) 

5 (OFF-FIL) 

6 (ON-FIL) 

7 (HL-A) 

8 (M12off) 

9 (LX-3) 

10 (NEP1) 

11 (NEP2) 

12 (LL21) 

13 (LL10) 
Rl M12on h 
R2 MP235 



2.4±o.6 

3.6±o.5 

5.6±o.5 

3.0±o.4 

5.0±o.5 

5.0±o.5 
1 q+0-4 

1 - J ~0.5 

8.8±o.5 

1+0.3 
-0.5 
4+0.8 
-0.7 

4.6±0.5 

2.6±o.e 



1.9 



-1.6 
-1.5 



3.1: 

3.4 H 



1.2+0.4 
9.1 + 1.3 
2.1+0.5 

4.9±2.i 
2 9+ 21 

7 q + 1.2 
' -°-1.3 

13.41?-? 

o n+2.6 
O.Z_2 7 

1.7±0.7 



11.8] 



1 : 

6.25~+S (4.0±g;S) 1 



14.9±5.4 

18.3^ 

2+3.2 
+45 

31.5_^ 7 

4 -'-3.7 
13.7±4.7 

24.3^-2 
ll.Qili 

1+3.5 

+73 

25.1 y 2 

27.3t 7 ^ 

10.5±5.6 

19.5^ 8 5 



0.568 H 



29.r 



U.ODU_ 01Q 
;+0.008 
'-0.006 

U.ODZ_ 002 

n 566+ 002 

U.OOO_g 004 

0.564+0.004 
0.564±o.oo4 
0.570±o.oo4 
0.567;";™? 
0.568l;fe 

-0.006 
-0.005 
+0.019 
-0.018 
-U.UUb 
-0.006 
+0.005 



0.568 
0.554 



°-"-i.o 
2.3+0.7 

2.5+0.7 

4.1+0.5 

9.5±o.7 
k o+i.o 

J - z -0.6 

4.5+0.6 
4.4+0.5 
6.0+0.7 
8.9±o.5 
7.0±i.i 

+0.8 
0.9 
+ 1.3 
-0.7 



6.4 
1.7 



u.uu-±_ gig 

fi^4+ 010 
U.D04_o.oi4 

0.654s 
0-659±8;gi5 

o-653±S:88t 

fifi0 +0 - 008 

U.DDU_o Q06 

u.di(_ on 
fifi2 +0 - 011 

U.DDZ_o.olO 

0.654+0.006 

+0.004 
0.002 
+0.006 
J -0.008 

0.652±o.oio 

+0.022 
-0.034 



0.653 
0.656 



0.642 



0.9±o.4 
<0.4 
0.5±o.3 
2.5±o.4 

1.6±0.3 
1.6±0.4 

0.6±o.3 

1.8±0.4 

o 0+0.3 
2.6±0.7 

1.8±0.5 



0.4 



+0.3 
-0.3 



80.06/80 

93.58/97 
124.53/ 92 
154.66/118 
113.18/113 
134.71/111 

79.05/85 
206.76/192 
165.23/116 

56.20/51 
108.14 /93 

80.63/87 



0.564 
0.578 



2.9 ±0.5 

9 i +0.8 



0.654S 



0.673 



-0.007 



<0.4 141.69/135 
0.9±o.3 77.18/83 



a The absorption column densities for the CXB and TAE components were fixed to the values tabulated in 
Table 3. 

The temperature of TAE and SWCX+LHB components were fixed to the best fit values tabulated in Table 3. 
The O abundance of the both components was set to in the fits. 
c The unit is phtons s _1 cm -2 keV _1 str _1 @lkeV. 

d The emission measure integrated over the line of sight, (l/47r) J n e nnds in the unit of 10 14 cm~ 5 str -1 . 
e The unit is LU =phtons s~ cm -2 str -1 . 

1 Contribution of Ovu emission in the Gaussian fitting function is subtracted using K Q intensity (see text). 
s The centroid energy was fixed to obtain the upper limit of the intensity. 

h An additional spectral component was included to represent the point source. See table note (g) of 
Table 3. 

1 "Narrow bump" component. The O abundance of this component was also set to 0. See table note (h) 
of Table 3. 



from the Gaussian normalization of Ovm K a . The ratio, 
u, of Ovu K/3 to Ovu K a intensities is a slow function 
of the plasma temperature for thermal emission. As we 
will show later, the temperature of the TAE component 
determined from O emission is in the range of kT = 0.19 
to 0.23 kcV. Then u = 0.056. If the emission is due to 
SWCX, u = 0.083 (Kharchenko et al. 2003). We estimated 
the Ovu Kp intensity assuming that 2 LU of OvnK a is 
from SWCX and the rest from thermal emission of kT = 
0.19 to 0.23 kcV (see below). The difference of the estima- 
tion from the totally SWCX or totally thermal cases is at 
most 0.2 LU, which is smaller than the statistical errors 
of Ovm intensity. 

In Table 5 we summarize the results of the spectral 
fits. The centroid energy of Ovm had to be fixed to the 
theoretical value for LH-2, because we obtained only an 
upper limit for this data set. Except for this case, the 
centroid energies are consistently determined to be equal 
to their expected values within the systematic errors of 
energy scale calibration (5 eV, Koyama et al. 2007, also 
see the Suzaku Technical Descricription). 

The Ovu emission intensities of NEP1 and NEP2 are 
not consistent with each other within the 90 % statistical 
errors, while the Ovm emission intensities are consistent. 



The difference between the best fit values of Ovu is 1.9 
LU. In the previous section, we discussed that time vari- 
able geocoronal SWCX of 1.5 - 2 LU level can remain in 
our spectra. Thus, the difference can be partly due to 
contamination of the geocoronal SWCX. Temporal varia- 
tions of Hcliosphcric SWCX could also contribute to the 
difference. 

3.2. Correlation between Ovu and Ovm intensities 

In Figure 3, we plot the derived Ovm emission intensity 
as a function of the Ovu intensity, and include the data 
points for M12on and MP235 from Masui et al. (2009). In 
the figure we notice two remarkable characteristics. First, 
there is a floor in the Ovu intensity at ~ 2 LU. Second, all 
the data points except for three (HL-B, MP235, and LL10) 
approximately follow the relation (Ovm intensity) = 0.5 
x [ (Ovu intensity) - 2 LU]. Masui et al. (2009) proposed 
that the Ovm emission line of the midplane field, MP235, 
is associated with a ~ 0.8 keV component that may arise 
from faint young dM stars in the Galactic disk. As we 
discuss later, the low latitude field at b — 10°, LL10, is 
also considered to contain emission of similar origin. Such 
high temperature emission produces little Ovu. We thus 
consider that the floor of Ovu emission and the strong 




Fig. 3. Relation between Ovil and Ovm surface brightnesses for the 14 sky fields observed with Suzaku. The horizontal and vertical 
bars of data points show the 1 a errors of the estimation. The contribution of Ovil emission in the Gaussian function for Ovm 
K Q is subtracted (see Table 5 and text). The diagonal lines show the relation between Ovil and Ovm, assuming an offset Ovn 
emission of 2.1 LU and emission from a hot plasma of the temperature and the absorption column density shown in the figure. The 
Galactic absorption column density of the observation fields are indicated by the maker size of the data points. The short names of 
five data points on the intensity floor of Ovil emission are also shown. 



correlation between additional Ovn and Ovm intensities 
suggest that the Ovn emission consists of two emission 
components of different origins; an emission component 
which emits Ovn of ~ 2 LU and no or weak Ovm emission 
(the floor component), and the other component which 
emits Ovn of ~ to ~ 7 LU and the Ovm intensity is 
about a half of Ovil (the linear component). 

Taking into account possible systematic error of Ovn in- 
tensity due geocoronal SWCX (~ 1.5 LU), M12on and LH- 
2 are also on the floor. The Ovn emission of the M12on 
field must arise from the near side of the MBM-12 molec- 
ular cloud whose distance is 60 - 300 pc. Assuming the 
midplane space-averaged neutral density (Ferriere 2001), 
the absorption length of Ovn emission is estimated to 
be 300 pc for MP235. The total Galactic absorption col- 
umn density ofLLlO (|6| = 10°) is 2.7 x 10 21 cm~ 2 , through 
which only 10 % of Ovn photons can penetrate. These 
indicate that the floor component is "local" emission and 
is associated with the SWCX + LHB component of the 
spectral model (see also discussion in the first paragraph 
of section 4.2). The emission temperature of the SWCX 
+ LHB component is kT ~ 0.1 keV when it is described 
with thermal emission. This is consistent with the weak 
Ovin emission of the floor component. 



Then the linear component is likely to arise from more 
distant part of our Galaxy. It is remarkable that LL21 
(|6| = 20°) contains significant linear component (6.4 - 2 
= 4.4 LU), even though the low Galactic latitude. The ab- 
sorption column density of this direction, 7.2 x 10 20 cm -2 , 
is about quarter of LL10 column density and the trans- 
mission of the Ovn emission is about 54 %. This sug- 
gests that a large fraction of the linear component arises 
from beyond the bulk of Galactic absorption, although the 
possibility that the difference between LL10 and LL21 is 
mainly due to spatial fluctuation of the emission measure 
cannot be ruled out. We thus consider that the linear 
component is associated with the TAE component. This 
component emits Ovn and Ovm with an intensity ratio 
of about 2 to 1 and shows field-to-field variations of ~ to 
~ 7 LU in Ovn intensity. The lines in Figure 3 show Ovn 
- Ovm intensity relations for different temperatures and 
absorption column densities. We have assumed 2.1 LU as 
the Ovn floor intensity, using the value for the midplane 
field (MP235). The lines suggest that the average tem- 
peratures along each line of sight of the plasmas emitting 
the Ovm emission and the Ovil emission above the floor 
are in a relatively narrow range of kT ~ 0.19 to 0.23 keV 
if they are arising from collisionally equilibrium plasma 
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in spite of the large Ovn intensity variations (~ to ~ 7 
LU). 

In Figure 4, we plotted the Ovn and Ovin emission in- 
tensities as functions of absolute Galactic latitude |6| (a), 
and as a function of absolute ecliptic latitude (b). The 
two data points marked with an open circle are M12on 
and MP235. The correlations of Ovn and Ovin emission 
intensities to Galactic latitude and ecliptic latitude are 
not very clear. Near solar minimum, Hcliosphcric SWCX 
emission is expected to be enhanced in low ecliptic latitude 
region of < 20° by ~ 1.5 LU (Koutroumpa et al. 2006). 
Such small difference may be masked by the spatial varia- 
tions of the TAE, systematic errors in geocoronal-SWCX 
removal, and time variations of Heliospheric SWCX. The 
relevance of Galactic latitude (for TAE) and ecliptic lat- 
itude (for SWCX emission ) are discussed in sections 4.2 
and 4.1, respectively. 

3.3. Spectral models consistently describe O emission 

If the Ovn emission originates in the two plasma emis- 
sion components of the model as discussed above, we 
would expect approximately constant temperature and 
normalization for the SWCX+LHB component that pro- 
vides the floor to the Ovii rate and a constant temperature 
(kT ~ 0.2 keV) for the TAE component with its constant 
Ovili/Ovn ratio. However, the results of the spectral fits 
shown in Table 3 (and Table 4 do not show such tendency. 
The Ovii emission intensities of the SWCX+LHB compo- 
nent estimated using the best-fit parameters in Table 3 
vary from 2 to 5.5 LU. The temperature of the TAE com- 
ponent is significantly higher than 0.2 keV and is far from 
constant. 

We consider that this discrepancy is related to the Fe 
and Ne emission. When we fixed the abundances (Fe to 
O, and Ne to O ratios) of the hot gas in the spectral fits, 
the temperature of the TAE component was mostly de- 
termined by Fe-L and/or Ne-K emission rather than O 
emission. Then the temperature and the intensity of the 
SWCX+LHB component was optimized to fill the remain- 
ing emission. We notice relatively large residual near Ne K 
emission energy, in some of the residual plots in Figure 2. 
We will now try fixing the temperature and the normal- 
ization of the SWCX+LHB component, and allowing the 
Ne and Fe abundances of the TAE component to vary to 
see if we can still get good fits and if the temperatures of 
the second component will be closer to 0.2 keV. This is 
model 2. 

As discussed in the previous subsection, the Ovn emis- 
sion from the three fields, M12on, MP235, and LL10 con- 
tains only the SWCX+LHB component. The temperature 
and normalization (emission measure) derived for that 
component (see Tables 3 and 4) are, respectively, about 
0.11 keV and 14 x 10 14 cm~ 5 str -1 for all three fields, which 
is consistent with the parameters obtained for model fit- 
ting of data from the LMC X-3 vicinity (LX-3) by Yao et 
al. (2009) (0.103 keV and 18.4 xl0 14 cm" 5 str" 1 ), which 
were derived from previous XMM-Newton and Suzaku 
observations of directions with absorption gas. For the 
model 2 spectral fits, we therefore adopted the best fit 



value of the midplanc observation (Table 3) as represen- 
tative of SWCX+LHB parameters. 

In Table 6 and Figure 5 we show the results of the fits. 
The x 2 values are generally good with the reduced \ 2 
values in the range of 0.95 - 1.28 for the degrees of freedom 
of 91 - 196, although most of them are larger than those 
for the fit with the first model in Table 3. For five of the 
spectra, the abundances were not constrained well. We 
thus fixed the Fe and Ne abundances to the solar value. 
For LH-2, the existence of the TAE was not significant. 
We estimated the upper limit of the intensity fixing the 
temperature to the average value of the other spectra. 
Excluding HL-B and LL10, the temperature of the TAE 
component are in the range of 0.18 to 0.24 keV with an 
average of 0.22 keV, which is significantly lower than the 
temperatures obtained by model 1. 

In model 2, the strong Fc-L and Ne K lines are explained 
by over-abundance of those elements. In particular, ~ 3x 
solar abundance was required for Ne for four fields. We 
consider that the strong Ne and Fe emissions could be 
also represented by a higher temperature emission com- 
ponent with solar abundances. For model 3, we fixed the 
temperature and the abundance of the TAE component 
to the average of model 2 (0.222 keV) and to the solar 
value, respectively. We introduce a fourth emission com- 
ponent with higher temperature which we denote TAE'. 
In Table 7, we show the results for the four fields. The 
resultant \ 2 values are comparable to those for model 2. 
In model 3, the excess Ne and Fe emissions are explained 
by emission with a temperature in the range of 0.6 to 0.9 
keV, and an emission measure of (1 — 2) x 10 14 cm~ 5 str -1 . 

4. Discussion 

4-1. Heliospheric Solar wind charge exchange 

Although we have been careful to avoid possible con- 
tamination by geocoronal SWCX through the selection 
procedures discussed in section 2, contributions from 
SWCX in interplanetary space are at best only partially 
removed by these methods since the long transit times 
through the solar system wash out the very large short- 
term intensity variations seen by monitor satellites located 
in the ecliptic plane near the Earth. Models of this he- 
liospheric SWCX emission are highly uncertain, but some 
predict that a large fraction of the Galactic R45 emission 
is from this source (Koutroumpa et al. 2006; Koutroumpa 
et al. 2008), making it one of the most severe limitations 
in determining the true interstellar and halo contributions. 

Solar activity affects the solar wind densities and ioniza- 
tion temperature. Near the solar maximum, the slow solar 
winds which have high ionization temperatures and high 
densities are ejected from the sun resulting in stronger 
SWCX induced Ovn and Ovin emissions (Koutroumpa 
et al. 2006) than in solar minimum. Near the solar mini- 
mum, slow solar winds are emitted from the equator region 
of the sun, and high speed, low-density, low-ionization- 
temperature winds are emitted from the high latitude 
region of the sun. Discrepancies between early XMM- 
Newton or Chandra observations which were also made 
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Table 6. Results of three-component (CXB, TAE, SWCX+LHB) spectral fit with SWCX+LHB temperature and normalization 
fixed, and with double broken power law CXB (model 2) 



ID 


iV H a 


CXB b 




TAE 






X 2 /dof 




10 20 cm- 2 


Norm c 


kT (kcV) 


Ne 


Fe 


Norm d 




1 (GB) 


1.40 




222 +u - in 


l c 


1° 


l- z Ho 


95.12/ 92 


2 (HL-B) 


3.36 


o c+0.5 
■^-J — O 6 


296+° : ° 93 

o.049 


1° 


1° 


1-618:? 


95.14/88 


3 (LH-2) 


0.56 


q c+0.4 
°-°-0.5 


0.222 


1° 


l c 


< 0.7 


95.33/100 


4 (LH-1) 


0.56 


5 6 +0 - 5 


fl OO7+0.262 
U.zo ( _o.058 


o 7 n+6.43 
°- ' y -2.02 


2 4 n+13.60 
z -^ y -1.98 


2.0l° ; 6 6 


162.31/127 


5 (OFF-FIL) 


4.19 


q 1+0.4 
J - 1 -0.4 


n 181 +0.008 
U - i8i -0.008 


o co+1.01 
°- do -0.80 


o qq+ 2 - 75 


15-311:1 


159.33/122 


6 (ON-FIL) 


4.61 


5 +0 - 5 

°- u -0.5 


U.zo»_ 2 5 


3.46ll;? 9 


1 ««+ i:97 
J-^-o.ss 


r 9 +0.9 
°- z -0.9 


130.94/117 


7 (HL-A) 


1.02 


^• y -0.5 


o 24 o+0.031 
U.z<±o_ 032 


0.91™ 


o en +1 - 16 

U - O1 -0.48 


4.318:1 


138.89/115 


8 (M12off) 


8.74 


2 O+0.4 
z - u -0.4 


U -184_ 030 


l c 


1° 


4.91M 


89.83/91 


9 (LX-3) 


4.67 


O n+0.5 

8-O-0.5 


n 91 O+0.026 

u - zlo -0.018 


3 f)4 +1 - 282 
°- u +-1.168 


o i 9+1.978 
°- lz -1.590 


7 7+1.2 
'■'-0.5 


210.94/196 


10 (NEP1) 


4.40 


q O+0.5 
J - z -0.4 


1Q1 +0-006 


2.09i°;f 7 


1 fiQ+0-75 


15 3+ - 9 

io -°-0.9 


164.60/120 


11 (NEP2) 


4.40 


o 4+0.7 


2 06+ 041 

U.ZUD_ gig 


2.65±i|| 


1 rq+1.743 

i.oy_i i 77 


10.811J 


57.94/55 


12 (LL21) 


7.24 


4 7 +o.5 

+ • ' -0.5 


iqq+o.019 

U.l»O_ 016 


2-ooioi^ 


z - 88 -1.33 


ii.Hl:? 


112.77/97 


13 (LL10) 


27.10 


2 c+0.7 


n 7 r 9 +0.240 


l c 


1° 


1-618:? 


91.44/91 



a The absorption column densities for the CXB and TAE components were fixed to the tabulated 
values. 

Two broken power-law model was adopted. The photon indexes below 1.2 keV were fixed to 
1.52 and 1.96 . The normalization of the former index component is fixed to 5.7, and only the 
normalization of the other component was allowed to vary. 



c The normalization of one of the broken power-law components with the unit of 
phtons s _1 cm~ 2 keV _1 str _1 @lkeV. 

d The emission measure integrated over the line of sight, (l/47r) / n c nnds in the unit of 
10 14 cm- 5 str- 1 . 

c The abundance was not constrained well, thus fixed to the solar value. 



Table 7. Results of four-component (CXB, TAE, TAE', SWCX+LHB) spectral fit with SWCX+LHB normalization and tempera- 
ture , and TAE temperature fixed, and with double broken power law CXB (model 3) 



ID 


10 20 cm- 2 


CXB b 
Norm c 


TAE 

kT (kcV) Norm d 


TAE' 
kT (keV) 


Norm d 


X 2 /dof 


4 (LH-1) 

5 (Off-FIL) 

6 (On-FIL) 
9 (LX-3) 


0.56 
1.90 
9.60 
4.67 


c q+0.5 
°- ,i -0.6 

2 6 +0 - 5 

z - D -0.5 
c 1+0.5 
°- 1 -0.6 

8 5+ - 7 

°-°-0.6 


0.222 
0.222 
0.222 
0.222 


i.9l u o: b 6 
10-318:? 

9.011-] 
' -^-i.i 


D 74R+0-326 

o.86ii8:ig 
o.676i8:8§S 

n ^^q+o i09 
u.ooy_ 193 


i-ol81 
o-8l8: 4 

i q+0-6 

1 «+°- 4 


160.48/128 
197.26/123 
122.93/118 
216.50/197 



a The absorption column densities for the CXB, TAE, and TAE' components were fixed to the 
tabulated values. 

Two broken power-law model was adopted. The photon indexes below 1.2 keV were fixed to 
1.52 and 1.96 . The normalization of the former index component is fixed to 5.7, and only the 
normalization of the other component was allowed to vary. 

c The normalization of one of the broken power-law components with the unit of 
phtons s _1 cm -2 keV _1 str _1 @lkeV. Add 5.7 to obtain the total flux at 1 keV. 

The emission measure integrated over the line of sight, i.e. (l/47r) J n c nnds in the unit of 
10 14 cm- 5 str" 1 . 
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Fig. 4. Ovn and Ovni emission intensities as functions of absolute Galactic latitude |6| (a), and as functions of absolute ecliptic 
latitude (b). The two data points marked with an open circle are Ovn intensities of the MBM 12 on cloud direction and the midplanc 
direction (235, 0). Both are taken from Masui et al. (2009). 



near the solar maximum in 2000-2002 and Suzaku observa- 
tions near the solar minimum were noted by Koutroumpa 
ct al. (2007) and Henley & Shelton (2008). For exam- 
ple, Lockman Hole (£,b) = (149.1,53.6) by XMM-Newton 
observations taken in October 2002 show much higher 
Ovn intensities (7 to 18 LU, Koutroumpa et al. (2007)) 
than those of Suzaku (2.5 and 4.1 from the present LH-1 
and LH-2 fields). Henley & Shelton (2008) showed that 
Ovn line intensities on and off directions of the shadow- 
ing filament (On-FIL and Off-FIL fields) were 10.65±g;^ 
and 13.86t};!g LU for thc XMM-Newton observations in 
2002 March, while they reported that thc intensities were 
6.51 and 10.53±g;|| LU for thc Suzaku observations 
(statistically consistent with the present results, 5.2lg'g 
and 9.5 ± 0.7 LU). If this difference is due to the different 
SWCX intensity as the authors suggest, it was brighter 
by about 5 LU during the XMM-Newton observations. 
The Ovn cmision intensity in the MBM12 on-cloud di- 
rection determined by Chandra (1.79 ±0.55) and Suzaku 
2.93 ± 0.45 were marginally consistent within the 90 % 
statistical errors, although the Ovni emission intensity by 
Chandra (2.34 ±0.36 LU) was larger than that of Suzaku 
(0.30 ±0.20 LU). 

The ROSAT sky survey was able to remove SWCX con- 
tributions with time variations on scales of a day or less, 
which should provide efficiencies in removing geocoronal 
SWCS similar to the procedures used here. The steadier 



part of thc hcliospheric contribution is still present, how- 
ever, and we might expect it to be considerably larger, 
since the ROSAT survey was conducted at solar maximum 
and the Suzaku observations near solar minimum. Both 
Suzaku and ROSAT are in low Earth orbit, and observing 
directions of the both satellites are almost perpendicular 
to the Sun-Earth-line. 

We calculated the R45 band counting rate expected 
for the best-fit model parameters of model 1 ( Table 3) 
and model 2 ( Table 6) . We used the ROSAT re- 
sponse function, "pspcc_gainl_256.rsp" , in CALDB at 
NASA/GSFC 5 . Two models gave consistent results within 
thc lcr statistical errors. 

Using the database at NASA/GSFC 6 , we extracted 
ROSAT R4-5 band average counting rate in a circular sky 
region of 36 arc minute diameter centered at the Suzaku 
XIS aim point. The size of the sky region is larger than 
than the Suzaku field of view. However, statistical pre- 
cision of survey will be inadequate on smaller fields. In 
Figure 6, we plot the observed ROSAT counting rate as 
a function of expected counting rate from thc Suzaku ob- 
servation. In the figure we include M12on and MP235 
fields using the results from Masui et al. (2009). The 
horizontal error bars in this figure is 1 a statistical errors 



available from ftp:/ /legacy.gsfc. nasa.gov/rosat/calib_data/pspc/cpf/matriceS/ 
available from http://heasarc.gsfc.nasa.gov/cgi- 

bin/Tools / xraybg/xraybg.pl 
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Fig. 5. Results of spectral fits with SWCX+LHB component fixed (model 2) except for MP235 (last panel). Observed spectra 
(crosses), best-fit model and its components (step functions) convolved with the instrument response function and residuals of the 
fit (bottom panels) are shown. The vertical error bars of data points correspond to the 1 a statistical errors. The model parameters 
of the SWCX+LHB component was taken from the midplane observation, MP235 in Table 4 
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Table 8. Comparison with ROSAT all sky survey 



ID 


Suzaku a 


RASS b R45 band 




10" 6 


cs -1 amin -2 


1 (GB) 


85.4 ±2.6 


186±21 c 


2 (HL-B) 


65.6 ±2.2 


150±21 c 


3 (LH-2) 


68.4 ±1.9 


85 ±15 


4 (LH-1) 


103.4 ±2.0 


103 ±16 


5 (Off-FIL) 


130.4 ±2.3 


227±34 c 


6 (On-FIL) 


105.3 ±2.1 


118±27 c 


7 (HL-A) 


103.7±2.2 


112 ±17 


8 (M12off) 


60.0 ±1.8 


99 ±20 


9 (LX-3) 


140.0 ±2.4 


180±16 c 


10 (NEP1) 


129.6 ±2.4 


158 ±4 


11 (NEP2) 


119.8 ±4.0 


158 ±4 


12 (LL21) 


109.3 ±2.8 


123 ±14 


13 (LL10) 


31.6 ±1.8 


44 ±9 


Rl (M12on) 


20.8 ±0.9 


37±11 


R2 (MP235) 


58.0 ±2.3 


109 ±21 



The errors in this table are all at 1-er signifi- 
cance. 

a The intensity expected in the ROSAT R45 
band from the Suzaku best-fit model function. 
Average intensity for a circular region of 

36' diamiter centered at the Suzaku-observation 
aim point. 

c These directions are considered to be contam- 
inated with the long-term enhancement: data 
points marked with open boxes in Figure 6. 

estimated from the spectral model fits of Suzaku data: 
we searched for the minimum and maximum expected 
ROSAT counting rate for the combinations of parameters 
on the y 2 = y„. + 1 surface in the 5 dimensional model 

"V- 'Wlllll 

parameter space. 

The observed and expected counting rates are well cor- 
related. However the observed ROSAT counting rates 
are systematically larger than the rates expected from 
the present Suzaku best-fit model. The sky directions 
of the five data points marked with an open rectangle 
(GB, HL-B, Off-FIL, On-FIL, and LX-3 ) are on streaks 
of bright areas which follow the RASS scan path. We 
thus consider that they are possibly contaminated with 
the long-term enhancement, i.e. with the SWCX from 
the geocorona, or by the solar X-ray scattering during 
the ROSAT observations (Snowden et al. 1994). If we 
exclude these data points, the scatter of the data points 
from a linear relation is significantly reduced. We then 
fitted the relation of the remaining nine data points with 
a linear function: (ROSAT observed rate) = p x (rate ex- 
pected from Suzaku) + q. We obtained ;? = 1.08±0.12 and 
q = (16.lii3'J) x c s_1 amm ~ 2 - The slope p is con- 
sistent with unity within the 90% statistical errors, while 
the offset q is positive. 

We consider that a large fraction of the positive offset 
can be explained by the difference in point source sensitiv- 
ity of the two sets of observations. For the ROSAT map 
used in the present analysis, point sources with count- 




Q 1^ I I I I I I I I I I I I I 

100 200 

Suzaku R45 intensity (10" counts s 'arcmin 2 ) 



Fig. 6. Observed ROSAT R45 band counting rate v.s. 
counting rate expected from Suzaku best-fit model functions. 
Points with open boxes are from possibly contaminated re- 
gions of the ROSAT survey. The vertical error bars are the 
1 a statistical errors of the ROSAT observation. The hor- 
izontal error bars are 1 a statistical errors estimated from 
the spectral model fits of Suzaku data. The diagonal solid 
line shows the predicted relation including the average point 
source removal correction, i.e. 12 X 10~ B c s — 1 amin -2 above 
the 1:1 relation. 

ing rate > 0.020 c s -1 in the R45 band were removed 
(Snowden et al. 1997). This threshold corresponds to 
3.8 x 10~ 13 crg cm -2 s -1 in 0.47 - 1.21 kcV assuming a 
power-law spectrum of a photon index of 1.95 and an ab- 
sorption column density of 5.6 x 10 19 cm -2 (see below). 
On the other hand, from the energy fluxes and their sta- 
tistical errors of the point sources we removed from the 
present Suzaku analysis, we estimate the typical detection 
threshold for point sources to be 1 x 10 _14 erg cm -2 s" 1 
in 0.47 - 1.21 keV. Thus there is a factor of ~ 40 differ- 
ence in the detection threshold. We estimated the aver- 
age surface brightness of point sources in the flux range 
of (1 - 38) x 10 -14 erg cm -2 s -1 in 0.47 - 1.21 keV, as- 
suming the log N — \ogS relations in 0.5 - 2 keV band 
obtained by ROSAT (Hasinger et al. 1996). To con- 
vert energy flux in 0.5 - 2 keV to the R45 band, we as- 
sumed a broken power-law spectrum with photon indices 
of 1.96 and 1.4 respectively below and above 1.2 kcV and 
various absorption column densities from 5.6 x 10 19 cm~ 2 
(Lockman hole field) to 1 x 10 21 cm~ 2 . The result we ob- 
tained was (14 — 10) x 10 -6 c s -1 amin -2 in the R45 band 
for N H = (0.056 - 1) x 10 21 cm- 2 . Thus the offset q is con- 
sistent with zero if we correct for contribution of point 
sources in the ROSAT data. Using the average value 
over the Nh range, 12 x 10 -6 c s _1 amin~ 2 , wc obtain 
q = (i.ltlli) x 10_6 c s_1 amin -2 . 

Although the ROSAT all sky survey was carried out 
near the solar maximum in 1990, the present results show 
that Suzaku and ROSAT i?^5-band intensities in at least 
nine directions were consistent with each other. The in- 
tensity of the spectral component for SWCX±LHB in 
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model 2, which contains Ovn emission of 2.1 LU produces 
a ROSAT R45 counting rate of 14 x 10~ 6 c s" 1 arnin -2 . 
Since this is comparable to the 90 % confidence upper limit 
of the offset q (= 17 x 10~ 6 c s" 1 amin~ 2 ). it suggests that 
on average, the maximum increase in hcliospheric Ovn 
emission intensity between the ROSAT and Suzaku is at 
most ~ 2 LU, or a factor of two. Since the solar maxi- 
mum increase should include lines from higher ionizations 
states such as Ovni and Neix that are more efficient in 
producing R45 counts, the increase in Ovn is probably 
much less than this. 

In section 2.3, we tried to remove as much as possible 
the time intervals in which the X-ray spectrum was con- 
taminated by the SWCX from the geocorona. Consistency 
between the ROSAT and Suzaku data also indicates suc- 
cessful removal. 

4-2. Intensity variations and temperature of the TAE 

Although the solar wind has a broadly simple struc- 
ture, with the slow wind confined to low ecliptic latitudes 
during solar minimum, the spatially non-uniform distri- 
bution of neutral H and He around the Sun leads to a 
complex pattern of SWCX emissivity. The distribution 
on the sky of observed SWCX emission is also a function 
of the observer's particular viewing geometry, as noted for 
example by Lallcment (2004) and Robertson & Cravens 
(2003), requiring detailed modeling which is beyond the 
scope of this paper. The general level of SWCX emission 
in the R45 band, however, is expected to be too small to 
explain the spatial variations observed in our work, con- 
sistent with our interpretation in model 2 that the spa- 
tially dependent components of Ovn emission are asso- 
ciated primarily with the TAE component and not with 
SWCX. 

Yao et al. (2009) constructed a thick hot disk model ex- 
tending above the Galactic disk in order to simultaneously 
explain the absorption and emission lines observed in the 
energy spectra of LMC X-3, the X-ray binary in the LMC, 
with Chandra, and in the energy spectra of the X-ray dif- 
fuse emission in the two directions about 30' away from 
LMC X-3 observed with Suzaku. They assumed density 
and temperature distributions exponentially decreasing in 
the direction perpendicular to the Galactic disk. They 
obtained as the best-fit parameters the scale heights of 
h T £ = 1.4(0.2,5.2) kpc and h n ^ = 2.8(1.0,6.4) kpc for tem- 
perature and density respectively, and the midplane gas 
temperature and H ion density of T = 3.6(2.9,4.7) x 10 6 K 
and no = 1.4(0.3,3.4) x 10~ 3 cm~ 3 , where £ is the filling fac- 
tor of the hot gas. The emission-measure weighted average 
temperature along a line of sight is T /(l + /i n /(2/ix)) ~ 
0.16 keV. This is lower than the best-fit temperatures of 
the TAE component with model 2 Table 6, although the 
discrepancy is not very large taking the the errors in To 
and h n /liT into account. 

For a plane-parallel configuration, we expect the in- 
tensity of the emission to increase from high to low lat- 
itude as oc sin -1 |6|,with a rapid decrease at low lati- 
tudes, b < 10°, due to Galactic absorption with a much 
smaller scale height than the emission. In Figure 7, in- 




20 40 60 

IGalactic latitudel (deg) 

Fig. 7. The emission measure of the TAE component mul- 
tiplied by sin | b| as a function of |6|. The TAE normalization 
factors are taken from model 2 (Table 6). The data point 
LL10 is not plotted because the higher temperature compo- 
nent of this field is like to have a different origin from other 
fields. 

stead of Ovn emission intensity, we show the emission 
measure of the TAE component for model 2 multiplied by 
sin | | (EM sin|6|). EM sin \b\ shows a large direction-to- 
direction fluctuation from a constant, which suggests that 
the hot gas is patchy and consists of number of blobs. 
The short angular scale spatial variation between the two 
Lockmann hole fields, LH-1 and LH-2, is puzzling. These 
are separated by only 0.42°, but the TAE component and 
Ovn emission are significantly stronger for LH-1 (Table 7, 
5). The ROSAT map also show a difference in R45 inten- 
sities of the two fields (Table 8), although the field of view 
of the two ROSAT fields are overlapped. Since the dis- 
tance between two lines of sight is only 10 pc at 1 kpc away 
from the sun, the density contrast of hot plasma must be 
high between outside and inside the blobs. Figure 7 may 
also indicate that EMsin\b\ is systematically larger for 
|6| < 50°, than for |6| > 50°. However, this could be a 
chance effect, and we need more data points, in particular 
for |&| > 50°. 

We estimate the parameters of hot gas by assuming an 
isotropic temperature of kT — 0.222 keV for simplicity. 
The total luminosity of the emission is estimated from 
the average value of EMsm\b\, assuming a plane parallel 
density distribution as 



L = 8n^A(T)EMsm\b\R z , 

where A(T) and R are, respectively, the emissivity per 
emission measure at a temperature T, and the outer cylin- 
drical radius of the emission region. Using the values of 
A(T) for bolometric flux (Sutherland & Dopita 1993) and 
0.3-2 keV band (APEC), we obtain, 



L b oi = 3.8 x 10 39 erg s" 1 



EMsin\b\ 



L 



0.3-2kcV 



3.6 x 10 14 cm" 5 str" 1 



1.1 x 10 39 erg s" 1 



R 



15 kpc 
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EMsia\b\ 



3.6 x 10 14 cm" 5 str" 
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The midplanc density and the total mass of the hot 
gas, Mtotj can be estimated by further assuming the scale 
height. 



n = 1.3 x 1CT 3 cm" 3 



EM sin I b\ 



1/2 



3.6 x 10 14 cm- 5 str" 1 



1.4kpc 



-1/2 




EMsm\b\ 



1/2 



3.6 x 10 14 cm- 5 str" 1 

\ 1/2 



h4 



1.4kpc 



As the scale height, h, of hot gas we assumed the tem- 
perature scale height from Yao et al. (2009), since it is 
smaller than the density scale height. 

Yao et al. (2009) concluded that the plasma responsible 
for the Ovil and Ovni absorption and emission towards 
LMC X-3 cannot be isothermal. However, the temper- 
atures averaged over line of sight for various directions 
determined by Ovil to Ovin ratio was remarkably con- 
stant (Figure 3). A possible origin of the hot gas is super- 
novae. In fact both the total mass and the total thermal 
energy of the hot plasma can be supplied at a supernova 
rate of 10~ 2 y _1 within the radiative cooling time of 3 
Gy. However, there is no particular reason why the gas 
should prefer kT — 0.2 keV. The spectra of the halo emis- 
sion of nearby star- forming galaxies and of some of normal 
galaxies are described by two-temperature models with 
fcT's in the range of 0.1 to 0.8 keV (Strickland et al. 
2004; Tiillman et al. 2006; Yamasaki et al. 2009). Thus 
we consider that kT = 0.2 keV is specific for our Galaxy. 
Then it is likely related to the virial temperature (kT = 0.2 
keV) for the rotation velocity of 200 km s -1 (e.g. Keres et 
al. (2005)). A possible explanation for this coincidence 
is that the hot gas was formed by cosmological accretion 
(Tfot et al. 2002; Rasmussen et al. 2009). The total 
luminosity of our Galaxy estimated above is consistent 
with the simulations. However, the radiative cooling time 
of the hot gas is 3Gy, and the hot gas needs be supplied 
at least on this time scale. Another interpretation is also 
possible that higher temperature gas has escaped from the 
potential and that the gas at the temperature near the 
virial temperature remains, since the escape time scale is 
by a factor of five shorter than the radiative cooling time 
scale. 

From the spectral fits with model 2, we obtained 2 to 
3 solar for [Ne/O] and [Fe/O] abundance ratios for four 
of the spectra. However, in model 3 (Table 7), we showed 
that these Fe-L and Ne " excess" emission can be explained 
by a higher temperature components of kT = 0.5 — 0.9 
keV. The HL-B field shows a high O emission temperature 
(0.30 keV, see Figure 3 and Table 6). Since all these fields 



0.1 



> 

u 
M 

bO.Ol 



io- 3 

2 

* o 



i' J i -- 
!_. i" 1 " 

'■■ : .J..^r' 
,.r l.j- M~ 

— L -H 1 KH — 1 — 


! '!-, 


"It 









0.5 



1 2 

Energy (keV) 



Fig. 8. Fit result of the Low Latitude 97+10 (data ID 13) 
spectrum using the faint dM star model by Masui et al. 
(2009) instead of the TAE component. 

are in high ecliptic latitudes, where only low-ionization- 
temperature winds are emitted from the sun during solar 
minimum, the Heliosphcric SWCX is not likely origin of 
the excess Ne and Fe emission. These lines of sight may 
contain blobs of high (> 0.22 keV) temperature hot gas 
which may be on the way to escape from our Galaxy. 

4-. 3. The spectrum of the b = 10° sample 

We obtained significantly higher temperature (0.75 
keV) for the TAE component for Low Latitude 97+10 
(LL10) than that of other samples (average = 0.2 keV). 
This direction has a high Galactic absorption of 2.78 x 
10 21 cm -2 , and the transmissions for Ovil and Ovin are 
respectively only about 10 and 20 %. Thus the emission 
from the thick disk of temperature ~ 0.2 keV will be sig- 
nificantly absorbed and hard to detect, if it exists in this 
direction. Masui ct al. (2009) detected an emission com- 
ponent of the similar temperature in the energy spectrum 
of the midplane direction, MP235. They suggested that 
the component is a sum of emission from unresolved faint 
dM stars existing between the bulk of Galactic absorp- 
tion and the Earth, and constructed model spectra assum- 
ing an average dM star spectrum with two-temperature 
thermal emission, the stellar X-ray luminosity distribu- 
tion functions, and the spatial densities of dM stars in 
the literature. The model spectrum could consistently ex- 
plain the observed spectrum not only in spectral shape 
but also in absolute intensity within 30 %. We fitted the 
LL10 spectrum using their dM star model spectrum con- 
structed for (£,b) = (90°, 10°), instead of the TAE com- 
ponent in model 1. As shown in Figure 8, this model 
reproduces the observed spectrum well (\ 2 = 89.2 for 90 
degrees of freedom). However, it was necessary to increase 
the intensity of the emission by a factor of about 5 from 
the model. This suggests that the emission from dM stars 
does not decrease so rapidly with increasing b as the model 
predicts, or that the there are large spatial fluctuations at 
b ~ 10°. We need more observations at low Galactic lati- 
tude in order to solve this problem. 
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5. Summary 

We presented spectra of soft diffuse X-ray emission in 
twelve fields observed with Suzaku together with the spec- 
tra of two other fields analyzed by Masui et al. (2009). 
In the data reduction, we carefully removed the contribu- 
tions of the solar-wind charge-exchange (SWCX) induced 
X-ray emission from the geocorona. However, ~ 1.5 LU 
uncertainty remains in Ovu intensity. The discrepancy 
between the two NEP observations (1.9 LU) can be partly 
due to yet incomplete removal of geocoronal SWCX. 

To determine O emission intensities, we first fitted the 
spectra with a model consisting of a broken power-law 
component for the CXB, and two thin thermal emission 
components of the solar abundance, one of which are sub- 
ject to Galactic absorption (model 1). The Ovu and Ovin 
emission intensities were determined using the model by 
setting O abundances of two thermal emission compo- 
nents to zero and adding two Gaussian functions to repre- 
sent Ovu and Ovin . The Ovu and Ovin intensities are 
strongly correlated and suggest the existence of an inten- 
sity floor for Ovu emission at ~ 2 photons s _1 cm -2 str -1 
(LU). The Ovin emission intensity of nine high-latitude 
fields show a tight correlation with the excess of the Ovu 
intensity above the floor. The relation is approximated 
as (Ovu intensity) = 0.5 x [(Ovu intensity) - 2 LU)]. 
These suggest that the Ovu emission arises from two ori- 
gins: approximately uniform emission of about 2 LU, and 
spatially variable (0-7 LU) emission from hot plasma of a 
temperature of ~ 2 keV. The former is likely to arise from 
Heliospheric Solar Wind Charge Exchange plus the local 
hot bubble (SWCX+LHB), and the latter from hot gas 
in more distant parts of the Galaxy, i.e. the transabsorp- 
tion emission (TAE). It is remarkable that for most of the 
fields, TAE average emission temperatures are confined in 
a narrow range (~ ±0.2 keV) around 0.2 keV. This tem- 
perature may be related to the virial temperature of the 
Galaxy, which locally corresponds to the rotation velocity 
of - 200 km s" 1 . 

The O emission intensities of the two thermal compo- 
nents of model 1 do not reproduce the above characteris- 
tics of O emission. This is because the abundance is fixed 
to the solar, while strong Ne and Fe-L emissions exist in 
some of the spectra. The observed spectra can be fitted 
with a model in which the intensity and temperature of 
the non-absorbed SWCX+LHB component are fixed to 
nominal values if we set the Ne and Fe abundances of ab- 
sorbed TAE component free (model 2), or if we include an 
additional higher temperature component (model 3). We 
found that four spectra required Ne to O abundances as 
large as 3 solar with model 2. Alternatively, these spectra 
can be fitted by model 3 with a higher temperature (0.5 - 
0.9 keV) emission component with solar abundances. The 
temperatures of the TAE component obtained with these 
two models were consistent with the values expected from 
the Ovu to Ovin ratio. 

The surface brightnesses estimated from the present 
best-fit model function were statistically consistent with 
the ROSAT All Sky Survey (RASS) map, even though 



the present observations were performed during solar min- 
imum, while the RASS was in solar maximum. The up- 
per limit for the difference was estimated to be 17 x 
10 -6 c s _1 amin _ i?^5-band counting rate after corrected 
for the difference in point source removal threshold be- 
tween the present Suzaku observations and the RASS. 

The origin of the TAE component was discussed in the 
context of the thick hot disk constructed by Yao et al. 
(2009). The emission measure determined by model 2 
shows a large deviation from the sin -1 |6| dependence ex- 
pected from a plane parallel configuration, which we con- 
sider to suggest short spatial scale structure of the hot gas. 
The total luminosity, midplane density, and the total mass 
of the hot gas were estimated assuming simple cylindri- 
cal geometry and a uniform temperature. The canonical 
temperature of the TAE component, kT ~ 0.2 keV, may 
be related to the virial temperature of our Galaxy. 

The lowest latitude sample of the present analysis, 
b = 10°, contains emission with kT = 0.75 keV instead 
of 0.2 keV. The spectral shape of this component can be 
represented by the faint dM star model by Masui et al. 
(2009). However the intensity must be increased by a fac- 
tor of about five from the model. 

The authors arc grateful to the anonymous referee 
for his comprehensive comments and useful suggestions, 
which improved the paper very much. They also would 
like to thank the Suzaku team for their effort for the oper- 
ation of the spacecraft, the calibration of the instruments, 
and the data processing. 
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